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Abstract. Helioseismology has shown that the chemical composition of the Sun has changed 
over its lifetime. The surface abundance of helium and heavy elements is believed to have 
decreased by up to 10% relative to their initial values. However, this reduction is too 
small to be tested by direct observations of the photospheric chemical composition. Here, 
we compare the predicted variations in the solar photospheric composition with precise 
measurements of abundances in meteorites and the solar wind composition. Although 
elemental composition ratios can vary by roughly a percent (e. g. for Ca/Mg and Ca/Fe) 
over the Sun's lifetime, their measurements are rife with uncertainties related to uncer- 
tainties in the interpretation of meteoritic measurements, photospheric determinations, 
and the complex fractionation processes occurring between the upper photosphere and 
lower chromosphere and the corona. On the other hand, isotopic ratios can be measured 
much more accurately and are not expected to be affected as much by extrasolar pro- 
cesses, although more work is required to quantify their effect. As the isotopic ratios evolve 
in the Sun proportionally to the mass ratios of the isotopes, light elements yield the high- 
est variations in isotopic ratios. They are predicted to reach as high as 0.6% for ^®0/^®0 
and are only slightly lower in the cases of ^^Mg/^^Mg and ^°Si/^^Si. Such a value should 
be well within the sensitivity of new missions such as Genesis. 



1. Introduction 

Modern solar models can be tested to a high degree of 
accuracy with newly available helioseismic data. Such mod- 
els infer the speed of sound in the solar interior by com- 
paring measured oscillation modes with model calculations. 
Because the speed of sound depends on the mean molecular 
weight of the gas (among other quantities), determining pro- 
files of sound speed yields information about possible com- 
positional inhomogeneities in the interior of the Sun. Indeed, 
solar models which include effects of elemental segregation 
succeed much better than those without in reproducing in- 
ferred speed profiles [Christens en- Dalsgaard et ai, 1996]. In 
these models, the photospheric abundances of all elements 
heavier than hydrogen have decreased by up to about 10% 
over the history of the Sun. However, this settling of the 
heavier elements out of the outer convection zone toward 
the center of the Sun has not been detectable with photo- 
spheric abundance measurements. In this work we consider 
the possibility of directly measuring this effect by measuring 
the composition of the solar wind. 

Lunar soils have preserved an archive of solar wind data 
dating back to the time of the formation of the lunar re- 
golith. Dust grains that have been exposed to the solar 
wind retain it with variable fidelity, although the lunar re- 
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golith as a whole appears to have preserved the entire so- 
lar wind fluence incident on the lunar surface [e. g. Wieler, 
1998]. Lunar soils are especially well suited to measure the 
isotopic and, to a limited extent, the elemental composi- 
tion of the noble gases and of few other elements, notably 
nitrogen. Early measurements of the He isotopic composi- 
tion appeared to show a secular variation of the solar wind 
''He/^He ratio [Geiss, 1973]. This could be interpreted as 
due to differences in the solar wind acceleration in the past 
or due to secular admixture of "^He from the "^He bulge in the 
Sun by a slow mixing process beneath the outer convection 
zone [Bochsler et ai, 1990]. The apparent secular change of 
the isotope ratio '^He/'^He has recently been re-interpreted. 
Heher [2002] analyzed lunar soil material of various antiq- 
uities for their noble gas isotopic composition. Apart from 
the known changes in He isotope composition, they found 
a change in the Ne isotope composition which is consider- 
ably larger than what could be attributed to the chemical 
evolution of photospheric composition. According to that 
work, the apparent secular change is due to slow removal 
of surface layers of lunar soil grains. The older grains have 
had more surface material removed, increasing the impor- 
tance of the deeper-seated and isotopically heavier so-called 
"SEP" or "HEP" component [Wieler et ai, 1986; Wimmer- 
Schweingruber and Bochsler, 2001a]. Heber [2002] used the 
measured correlation between Ne and He isotope composi- 
tion to account for this removal by assuming that the Ne 
isotope composition did not change over time. Then the 
apparent secular variation of the He isotopes vanishes and 
hence the variations of '^He/'^He are not necessarily a test 
for solar evolution models. Nitrogen is known to have large 
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secular variations in its isotopic composition in lunar soils 
[see e. g. Kerridge, 1993, for a review.]. However, recent in- 
vestigations by Wieler et al. [1999] have shown that about 
90% of the implanted nitrogen is non solar, i. e. due to some 
other, as yet unidentified, source. 

Because there is no unambiguous evidence in lunar soils 
for secular changes in solar wind composition which could 
be interpreted as due to elemental segregation, another ap- 
proach needs to be chosen. Here we investigate the possibil- 
ity of comparing the composition of certain, well determined 
meteoritic abundance ratios with the same abundance ratios 
measured in the solar wind. The formation of the solar sys- 
tem about 4.6 billion years ago left compositional signatures 
in its constituent bodies and in different samples of these 
bodies. Compositional differences can be used as tracers 
for formation and alteration processes active mainly in the 
early solar system, but partially continuing to the present 
day. For example, volatile elements are extremely depleted 
in terrestrial samples and in most meteorites. Elements dif- 
fer in their chemical properties, resulting in variations of 
elemental abundance ratios in different separates of vari- 
ous solar-system samples. Because the Sun contains 99.9 % 
of the mass in the solar system, its composition serves as 
the reference for solar system composition and often solar 
and solar-system composition are used synonymously. This 
information comes from spectroscopic abundance determi- 
nations of the solar photosphere or from in-situ measure- 
ments of the solar wind. However, much of the information 
available today comes from studies of the composition of 
meteorites and of meteorite mineral separates or samples 
from other solar-system bodies. Based on the comparison 
of such investigations the composition of the bulk solar sys- 
tem has become known to a remarkable degree of accuracy. 
In fact, because of the close match in elemental abundances 
between solar and Cl-meteoritic values, CI meteorites are 
considered the most primitive material in the solar system 
besides comets or the Sun. Figure 1 illustrates this good 
agreement using data from the compilation of Grevesse and 
Sauval [1998]. The right-hand panel is a plot of the loga- 
rithm of the Cl-meteoritic abundance ratios of element X 
with respect to Mg, X/Mg, divided by its corresponding 
photospheric value plotted versus 50% condensation tem- 
perature [Wasson, 1985]. This quantity is often used to 
order elemental abundance trends in different classes of me- 
teorites. The comparison of photospheric and Cl-meteoritic 
abundance ratios is not reliable for condensation tempera- 
tures below about 400 K. For all other elements, with the 
exception of Li, the agreement between the two different 
abundance measurements (Cl-meteoritic and photospheric) 
is considered remarkably good. This is illustrated in the left- 
hand panel of Figure 1. It shows a histogram of the abun- 
dance ratios, (X/Mg) /(X/Mg) photo. The histogram is sym- 
metric around zero, indicating that there is no systematic 
disagreement between the two sources of solar abundances. 
The standard deviation of the histogram is 0.084 dex, cor- 
responding to a width of 21.5%. A Gaussian fit with that 
value for a is plotted for comparison. The fit is satisfac- 
tory, and the result does not change when using a different 
bin size for the histogram, e. g. the one shown in thin lines. 
Assuming the width of the histogram to be representative 
of the quadratic sum of the uncertainties of Cl-meteoritic 
and photospheric we subtract quadratically the maximum 
uncertainty of Cl-meteoritic measurements, 10%, to obtain 
the expected uncertainty of photospheric abundance deter- 
minations, ~ 20%. This value is less than that given by Del 
Zanna et al. (2001), 25% - 30%, who consider this to be the 



uncertainties in atomic properties alone. Moreover, uncer- 
tainties in the treatment of the thermodynamics of the line- 
emitting regions in the photosphere and lower chromosphere 
can be substantial [Holweger, 2001]. Nevertheless, some au- 
thors do give smaller uncertainties for certain elements with 
well known atomic parameters. Meteoritic uncertainties are 
generally smaller, on the order of 3% - 10% and are mainly 
due to possible sampling biases and will be considered later. 
Thus, from this discussion, we conclude that abundance ra- 
tios of refractory, i. e. non- volatile, elements offer the best 
choice of measuring the effect of elemental segregation. 

We consider the physics of elemental segregation in sec- 
tion 2, apply this to solar evolution in section 3, where we 
also give predictions for the evolution of some elemental and 
isotopic abundance ratios. We consider possible fractiona- 
tion processes which would influence measurements in the 
solar wind and investigate the potential of various solar wind 
measurements in section 4. 

2. The Physics of Chemical Evolution in 
the Sun 

The basic structure of the Sun consists of a dynamically 
stable inner sphere which extends to a little more than 70% 
of the Sun's radius and an overlaying convective shell which 
extends to the photosphere. The turbulent motions in the 
convection zone are so rapid compared to other time scales 
that one can assume that it is well mixed, leading to a homo- 
geneous composition in its entirety. Compositional changes 
due to nuclear processing in the convection zone can be ex- 
cluded because the temperature is too low. Hence changes in 
the chemical composition have to be generated at the bound- 
aries of the convection zone, either in the dynamically stable 
region at its base or at the solar surface (through infall of 
material or mass loss). As any excess or deficit of an ele- 
ment is diluted in the entire convection zone, a variation of, 
say, 10% in the observed abundance entails that the mass 
of the element i in question added or subtracted from the 
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Figure 1. Comparison of CI - meteoritic and photospheric 
abundances normalized to Mg. The right-hand panel is a 
plot of the logarithm of the CI - meteoritic abundance ra- 
tio X/Mg divided by its corresponding photospheric value 
plotted versus 50% condensation temperature. The left- 
hand panel is a histogram of the abundance values for two 
different bin sizes, also shown are a fit to the histogram 
and the uncertainty of the currently best solar wind abun- 
dance determination of N/0. The uncertainties (and hence 
the width of the histogram) are nearly entirely due to un- 
certainties in photospheric abundance measurements. For 
instance, atomic parameters limit measurements to an un- 
certainty of about 20 - 30 %. Meteoritic uncertainties are 
smaller, on the order of 3 - 5% and mainly due to possible 
sampling biases. 
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convection zone must be Mi — 0.1 Xi Mcz , Xi being the 
original mass fraction of element i and Mcz being the mass 
of the convection zone (in the current Sun Mcz is roughly 
2% of the total mass of the Sun [Christensen-Dalsgaard et 
al, 1996]). For instance, if one terrestrial planet had mi- 
grated into the Sun in the early stages of the formation of 
the solar system, this would have resulted in an increase 
of solar metallicity by only about 1% using the present-day 
mass fraction of the convective zone. Note that this process 
will not alter the relative abundances of refractive elements 
among themselves. Moreover, in the early Sun, the convec- 
tion zone was deeper and hence more massive and a larger 
flux of matter into it would have been necessary for any 
observable effect to be expected. 

The current mass loss or accretion rates are much too low 
to have any measurable impact on photospheric abundances. 
Earlier in the Sun's life it is probable that mass loss or accre- 
tion rates wore higher. Because of the faster rotation of the 
Sun in the past, dynamo action was enhanced. This made 
more energy in the form of magnetic fields available to drive 
the solar wind. Simple calculations show that this could at 
most have doubled the solar wind fluence over the lifetime 
of the Moon [Wimmer-Schweingruber and Bochsler, 2001b], 
commensurate with measurements of solar wind noble gases 
in lunar soils [Geiss, 1973]. Extrapolating these results to 
earlier times does not result in a significant increase of mass 
loss. Moreover, the convection zone was also much more 
massive in the past. During the Sun's evolution on the pre 
main sequence when mass exchange would be expected to 
be greatest the Sun is thought to have been completely con- 
vective [Hayashi, 1961]. The hypothesis of large mass loss 
rates in the young Sun has been examined but is not fa- 
vored by lielioseismology [Morel et al., 1997]. The addition 
of material that has been processed by nuclear reactions in 
the dynamically stable core can also be excluded. Again, 
temperatures are high enough only in its innermost regions 
making any exchange at its boundary impossible (with the 
sole exception of lithium, see section 2.2). 

2.1. Migration and Separation Processes 

The evolution of abundances in the solar convection zone 
and in the outer radiative envelope just below it occurs 
chiefly as a result of the differential effects of gravity on 
the elements of different mass. In a simple two component 
gas, constituted of hydrogen and helium, the inward grav- 
itational pull is four times greater on helium than on hy- 
drogen. As a result, helium will gradually settle toward the 
Sun's center and the fractional abundance of hydrogen will 
rise at the surface. Differential radiation pressure induces 
a minor correction in the overall chemical evolution in the 
convection zone but is plays a dominant role in the evolution 
of abundances ratios. 

The rate at which gravitational settling occurs is deter- 
mined in a largo part by the "friction" within the gas, which 
is dominated by Coulomb interactions. Because friction in- 
creases with the overall ionization state of the plasma, the 
settling rate generally will decrease as temperature and den- 
sity rise. 

As the evolution of photospheric abundances is governed 
by the evolution of the abundances at the base of the con- 
vection zone, the deeper the convection zone, or the smaller 
the stable radiative interior, the smaller are the fluxes of 
elements in and out of the convection zone. This is the 



combined effect of enhanced interactions between particles 
leading to smaller diffusion velocities and the geometrical 
effect of reducing the area of the surface through which par- 
ticles can flow. Moreover, the dilution of abundance changes 
in the larger convection zone further increases the time scale 
for the rate of change of composition. This is important be- 
cause the depth of convection zone varied strongly as the 
Sun evolved [as shown in Bahcall et al., 2001]. 

The time scale for the evolution of surface abundances 
can be approximated as [Michaud et al, 1976] 
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where p is the density, r is the radius, and Vdis is the dif- 
fusion velocity all evaluated at the base of the convection 
zone (CZ). The timcscales for helium, silicon and iron with 
respect to time shown in Figure 2 illustrate the combined 
effect of the increases in Wdifi and the radius of the base of 
the convection zone and the reduction in Mcz- 

The gravitational force exerted on an ion species "i" is 
AiTUH g{r), where g{r) is the gravitational acceleration at 
distance r from the center of the Sun. Radiative pres- 
sure is included in the modeling of diffusion by chang- 
ing the effective gravity felt by different atomic species to 
A-i niH {g — 5rad;i), whcrc (7rad:i is the radiative acceleration. 
ffrad;i is vcry scnsitivc to the to features of the species spec- 
trum, i. e. to the atomic structure of the ion, on the ioniza- 
tion state of the ion and of other ions competing for the 
same photons. Completely ionized elements feel practically 
no radiation acceleration. 

The expression for g^ad-.i is taken from Richer et al. [1998] 
(please refer to this paper for a detailed discussion) as being 
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Figure 2. Evolution of the timescale for the variation of 
chemical composition in the convection zone with the age 
of the Sun, the values are normalized to those in the con- 
temporary Sun at log(age)=0.66. The timescale is element 
dependent, here helium (t(Hc), normalized to 39 Gyr)), sil- 
icon (t(Si), normalized to 51 Gyr) and iron (t(Fe), normal- 
ized to 51 Gyr) arc shown. Also shown are the evolution of 
the mass of the convection zone {Mcz, normalized at the 
current value of 0.023 Mq), and the radius of the base of 
the convection zone {Rcz, which is currently at 0.710 R© 
in this model). 
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In this expression Xi is the mass fraction of species i, /aux 
is the fraction of the radiative flux which will be transferred 
to species i (refer to Eq.l in Richer et al. [1998]). The 
Rosseland mean opacity kr is a measure of the opacity of 
the plasma and determines the total momentum of the ra- 
diation field transferred to the absorbing plasma. 

The inverse relationship between g-cad-.i and the species 
abundance has an interesting potential consequences for iso- 
topic fractionation. Different isotopes of an element can 
have large differences in abundance, which would favor the 
less abundant species. This could potentially create large 
isotopic anomalies. The condition for such an effect is that 
the spectral lines of the minor isotope be well separated from 
the spectral lines of the major isotope and that there would 
be no coincidence with strong lines from other elements so 
that /flux can become non negligible for the minor isotope. 
This is not the case in the Sun, not even for helium for 
which the effect would be the most noticeable. The temper- 
atures below the convection zone where separation by radia- 
tive pressure could occur arc large enough, T > 2 X 10® K, 
that the thermal width of the spectral lines is much larger 
than the line separation for different isotopes of a given ele- 
ment. Therefore, /hux will be small and will cancel out the 
effect of the 1/Xi factor. 

It follows that in all calculations of the evolution of iso- 
topic ratios in the present paper the .grad;i will be assumed 
to be the same for all isotopes of a given clcnicnt. The 
radiation pressure will only play a role in the evolution of 
elemental ratios. 

2.2. Mixing below the solar convection zone 

The first evidence for mixing outside of the convection 
zone came from the strikingly low photospheric lithium 
abundance. Compared to its meteoritic abundance it is de- 
pleted by a factor of 137 in the photosphere [Carlsson et al., 
1994]. Solar models without mixing cannot account for this 
large difference. Modern solar models include various forms 
of mixing, based on convective overshoot or rotationally- 
drivcn mechanisms, and can now account for the observed 
lithium depletion while preserving the observed beryllium 
abundance and ^He/''He ratio [See e.g.Brun et al., 2000; 
Richard et al., 1996]. 

Moreover, based on helioseismic results, it appears that 
diffusion alone is not sufficient to explain the radial profile 
of the speed of sound in the Sun. There are large discrepan- 
cies between older models and measured values, especially 
around the lower boundary of the convection zone. This can 
be and has been rectified by including a limited extension of 
mixing below the convection zone into the stable radiative 
interior [Christensen-Dalsgaard, 1996]. 

Both issues, lithium and sound speed, can be resolved si- 
multaneously by the same mixing mechanism [as in Brun et 
al, 2000] The mixing occurs in a narrow region at the base 
of the convection zone in which there is a largo gradient in 
rotational velocity, this region is called the tachoclino. Al- 
though several mixing processes in that region can match 
the constraints in the Sun, studios in other solar type stars 
[Piau and Turck-Chieze, 2002] suggest that turbulent mix- 
ing caused by the shear due to the rotational velocity gradi- 
ent [Spiegel and Zahn, 1992] is superior to other processes, 
the most often suggested being convective overshoot [e.g. 
Blocker et al., 1998]. 

3. The predicted evolution of the composition 
of the outer convection zone 

For this work, we have used the models calculated by 
Turcotte et al. [1998] which feature the most detailed treat- 



ment of microscopic diffusion in the Sun. The evolution 
of the abundance of the minor isotopes of elements heavier 
than carbon have been calculated explicitly. The Turcotte 
et al. [1998] models have also been used by Bochsler [2000] 
to investigate the effect of isotopic fractionation in the Sun 
on solar wind abundances. He estimated the fractionation 
based on isotopic mass differences and the published factors 
of depletion for the major isotope of each element. In the 
current paper the full effect of the solar evolution on the net 
isotopic fraction is taken into account and an approximate 
prescription for the effect of mixing beyond the convection 
zone is also used. The net effect on isotopic ratios is quite 
similar in Bochsler [2000] and here in the absence of mix- 
ing. While the Turcotte et al. [1998] models do not include 
the effect of mixing beneath the convection zone, we have 
approximated how mixing affects the evolution of the photo- 
spheric composition by reducing the changes in composition 
such as to reflect the results of Brun et al. [2000] who include 
mixing in the tachocline in addition to a standard treatment 
of diffusion. 

Because of the assumptions on initial composition inher- 
ent in the computation of solar models, the absolute values 
of the abundances predicted by solar evolution models can- 
not be used, only relative changes are meaningful (with the 
notable exception of the helioseismic measurements of the 
helium abundance). Here we assume an initial composition 
as in Grevesse and Noels [1993] calibrated to reproduce a 
contemporary value of Z/X of 0.245. This value of Z/X 
inferred from the observations is thought to be accurate to 
10% at best. 

The predicted surface composition at 4.6 Gyr is illus- 
trated in Fig. 3 where the difference between the compo- 
sition at 4.6 Gyr and the original composition is plotted in 
function of the atomic number. The features reflect mostly 
the effect of the differential radiative forces. These effects 




5 10 15 20 25 
Atomic number 



Figure 3. Relative change in the surface abundances pre- 
dicted by the detailed solar models with diffusion, the lower 
one without mixing and the higher one with an approximate 
affect of mixing at the base of the convection zone. 
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are relatively small however as the spread in the relative 
abundance variations spans 8.8% for Ar to 7.6% for Ca. The 
mean variation in the detailed model is around 8.5%. The 
model including mixing has a lower relative change, with an 
adopted value for all estimates done in this work of 71.8% 
as high as that of the detailed diffusion-only model. 

As a consequence of the evolution of the Sun itself, re- 
flected in the decreasing timescales shown in Figure 2, the 
rate of abundance changes has increased sharply as the Sun 
has aged. This is illustrated in Fig. 4 where the elemental 
abundance ratio of calcium to iron and the isotopic abun- 
dance ratio of ^^C and ^^C are shown to increase exponen- 
tially in the near past. One notices that the trends are 
the same for the elemental and isotopic ratios even though 
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Figure 4. The evolution of the calcium over iron and 
i2(-i^i3(-i abundance ratios in the solar convection zone over 
it's lifetime. The solid line and dotted lines show Ca/Fe 
with and without mixing below the convection zone respec- 
tively. The dashed and dot-dashed lines show ^^C/^'^C with 
and without mixing. 



radiative forces are irrelevant in the latter case. One con- 
sequence of this rapid recent evolution is that fossil records 
of the solar wind at different epochs (lunar soils for exam- 
ple) may possibly contain evidence of the variations of some 
abundance ratios. 

The predicted evolution of isotopic ratios for He, C, N, 
O, Ne, Mg, Si, and Ar is presented in Table 1. Recalling 
the discussion on diffusion processes (Section 3), one can 
infer that the onlv difference of note for isotoDes in the con- 
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Figure 5. Ca/Fe vs. Mg/Fe elemental mass abundance 
ratios for various studies of Cl-chondritic samples. The ab- 
breviations next to the symbols refer to the authors of the 
corresponding studies: P,S,&Z (1981), Palme, Suess, and 
Zeh, (1981); A&E (1982), Anders and Ebihara (1982); A&G 
(1982), Anders and Grevesse (1982); W&K (1988), Wasson 
and Kallemeyn (1988); L&F (1998), Lodders and Fegley, 
(1998); M&S (1995), McDonough and Sun, (1995); P&B 
(1993), Palme and Beer, (1993); MPI Mainz, Spettel et al., 
(1993). We have also plotted the standard deviation of MPI 
measurements of 8 samples as well as the uncertainty in 
the unrealistic mean. Single measurement errors are smaller 
than this uncertainty, the scatter reflects variability in dif- 
ferent mineral separates of CI chondrites Orgeuil and Alais. 
The diagonal of small box shows the effect of solar evolution 
on photospheric abundances. The absolute position of the 
box is unimportant, it describes the relative change of the 
two elemental abundance ratios. 



Table 1. Values of the relative variation in isotopic abundance 
ratios over the Sun's life {A(A/B) /(A/ B)o) in percent. The value 
without mixing is taken from the Turcotte et al. [1998] models 
and the values with mixing are the same results eorrected by 
the estimated effect of tachoeline mixing as discussed in Brun et 
al. [1999]. 
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Figure 6. Compilation of the silicon isotopic composi- 
tion of various sources in the solar system. See Wimmer- 
Schweingruber et al. [1998] for sources and detailed discus- 
sion. 
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text of diffusion is their mass. It is assumed licrc tfiat all 
characteristics (i. e. ionization state, radiative forces) for the 
major isotope were shared by the minor isotopes with the 
exception of mass. 

As a result of the large mass difference between its iso- 
topes, helium has the largest variation in its isotopic ratio. 
Although the variation of the isotopic ratio is by far the 
largest variation of all ratios either isotopic or elemental, 
the study of helium is compromised by the lack of accurate 
measurements of the isotopic ratio of helium at earlier times 
in the Sun's life, as has been touched upon in the introduc- 
tion. For other elements, the isotopic ratios vary much less 
but not much more so than elemental ratios as seen in Ta- 
ble 2. Not surprisingly, the largest effect, after helium, was 
found for ^*0/^''0 and ^^'At:/*°Ai which have similar iso- 
topic relative mass differences. As a general rule, as the 
mass of the element increases the effect decreases gradually 
as the relative mass difference between isotopes decreases. 
Nonetheless, significant effects can be seen when compar- 
ing the more massive isotopes of Mg and Si to their major 
isotopes. 

4. Selecting the best tracers of the solar 
chemical evolution 

The aim of this work is to identify elemental or isotopic 
abundance ratios that are especially susceptible to elemental 
migration effects, well determined in meteorites, and which 
can be precisely measured in the solar wind. In this sec- 
tion we consider possible biases and fractionation processes 
which need to be taken into account in both metcoritic and 
solar wind studies and then motivate out selection. 

4.1. Elemental Abundances From Meteorites 

Because meteorites can be analyzed in the laboratory, el- 
emental abundance measurements performed on meteorites 
or terrestrial samples reach the highest possible precision. 
Measurement uncertainties are vanishingly small even for 
the application considered in this work. In Figure 5 we 
have plotted the Ca/Fc clement abundance ratio versus the 
Mg/Fe ratio as derived by different authors (indicated next 
to the symbols and referenced in the caption). The actual 
measurement uncertainties are considerably smaller than 
the spread in all the data and also smaller than the spread 
in the data from the Mainz group, who publish abundance 
data for 8 samples. Their reported uncertainty in the mean 
is about 1% for Mg/Fe and larger for Ca/Fe. It is important 
to note that these uncertainties are not primarily due to the 
different authors of the studies, but rather to different (some 
subjective) choices of samples of metcoritic matter. More- 
over, there are only five known CI chondrites, of which one 
has been studied extensively (Orgeuil), one in some detail 
(Alais), and three are so small that obtaining samples from 
them must be considered a feat on its own. Hence the ob- 
served variations in samples in Figure 5 must be considered 
a lower limit on the variability of these elemental abundance 
ratios in different separates of the most faithful samples of 
the presolar nebula. For comparison, we have also plotted a 
box with a diagonal which shows the maximum evolution- 
ary change of the photospheric abundance ratios based solar 
evolution models without mixing. Obviously, even if solar 
wind measurements were possible with much better preci- 
sion than currently possible (for instance with the Genesis 



samples) , the effect of solar evolution can not be determined 
by considering elemental abundance ratios. Uncertainties in 
metcoritic abundances will render futile any such endeavor. 

4.2. Isotopic Composition From Meteorites 

It is well known from studies such as the one of Anders 
and Grevesse [1989] that there is little variation in the abun- 
dance ratios among the non-volatile elements in different 
materials of the solar system. This homogeneity in elemen- 
tal compositions of different solar system materials reflects 
the fact that the early condensates in the inner parts of the 
solar nebula did not undergo strong chemical fractionation. 
Isotopic fractionation - e. g. due to incomplete condensation 
- was limited to tiny effects [Humayun and Clayton, 1995b]. 
Consequently, it is not surprising that there is also not much 
variation in the isotopic abundance of the refractory ele- 
ments in the bulk of solar system samples. 

The small variation of solar system isotopic abundances 
in refractory elements, is illustrated in Figure 6 for the 
case of Si. The conventional 5 notation is used, ^'^''Si 
= 1000 X ((n"»Si/n2SSi),ampie/(n''Si/n2»Si),rf - 1) in per- 
mill. In this three-isotope diagram a variety of different 
solar system samples follow a straight line with slope 1/2 
which is typical for chemical isotope fractionation. The en- 
tire variability is confined to a narrow range of few permills 
per mciss unit and hence the uncertainty of the mean of the 
~ 200 samples is much smaller, about .15 permill. This 
clearly reflects the efficient isotopic homogenization of the 
solar nebula. It also lends confidence to the assumption that 
there cannot be rrmch difference between the bulk isotopic 
composition of Si and other refractory or weakly volatile el- 
ements in the Sun and the rest of the solar system samples, 
including terrestrial ones. It is for this reason that Si and 
Mg are often used as standards for isotopic fractionation in 
the solar wind (to within few permill per mass unit). 

4.3. Elemental Composition of the Solar Wind 

The composition of the solar corona is not the same as 
that of the photosphere and the differences allow us to char- 
acterize various fractionation processes. There arc several 
properties of an individual atom or ion that allow for frac- 
tionation. In the strong gravitational field of the Sun, frac- 
tionation according to mass would be extremely strong if 
there were no other forces acting on the particles. Sim- 
plistically assuming an isothermal atmosphere between the 
photosphere and the chromosphere, a barometric equation 
would tell us that we should not be observing any heavy 
elements in the corona or the solar wind at all. The expo- 
nential nature of this equation thoroughly depletes any mass 
but the lightest. Obviously, this is not the case, and other 
forces must be acting upon the ions and atoms already in 
the lowest layers of the solar atmosphere. Among the best 
known is the process that, in some still badly understood 
way, separates neutral and charged particles. Physically, it 
is probably based on a competition of two time scales, the 
first ionization time (FIT) and some time scale that charac- 
terizes the removal or supply of neutral material. Because 
FIT is a model dependent quantity, it is more conventional 
to speak of this important fractionation effect as the FIP 
(first ionization potential) effect. The FIP is a well deter- 
mined atomic property for all important elements in the 
solar wind. 

In order to describe the FIP effect, it has become cus- 
tomary to use a double ratio of abundances. One normally 
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determines the solar wind abundance ratio of element X with 
respect to oxygen, X/0, and normalizes this observed quan- 
tity with the corresponding photospheric value, X/Ophoto- 
The double ratio H = (X/0)/(X/Ophoto) lies between 2 and 
3 for low-FIP elements (FIP < 10 V) in the slow wind and 
between 1 and 2 in the fast wind [e.g. von Steiger et at, 
2000]. Because volatility is generally correlated with FIP 
this means that refractory elements (i. e. low-FIP) are en- 
hanced by a factor of about two in the average solar wind 
when compared with their photospheric abundance. 

Because of the magnitude of the FIP effect, only high- 
speed solar wind should be used for testing solar evolution 
models. However, even in this otherwise benign type of solar 
wind, elemental abundances are not necessarily fully deter- 
mined by the FIP effect. Especially wave-particle interac- 
tion, widely believed to be the accelerating agent in this type 
of solar wind, is highly dependent on the ionization state of 
the ion and hence on its atomic structure. Given current 
uncertainties in atomic parameters, it appears uulikoly that 
elemental abundances will be used in the near future for the 
purpose studied in this work. 

4.4. Isotopic Composition of the Solar Wind 

The strong fractionation observed for elements in the slow 
solar wind is mostly due to the FIP effect (First Ionization 
Potential). Although this effect is of utmost importance for 
elemental fractionation, it is probably quite unimportant for 
isotopic fractionation. For instance, the elemental fraction- 
ation model of Marsch et at. [1995] can be understood as 
a competition between two time scales, diffusion and ion- 
ization. The authors derive the following expression for the 
mass-dependent part of elemental fractionation fij between 
two species i and j. 



Aj + l A.J 
A-j + l 



1/4 



where Ai^j is the atomic mass number of species i and j. For 
the example of magnesium this expression evaluates to a de- 
pletion of the heavy isotopes of less than half a permill per 
meiss unit. In other words, FIP is unimportant for isotopic 
fractionation according to the Marsch et al. [1995] model. 
Note however, that the model of Schwadron et al, [1999], 
which invokes resonant ion-cyclotron heating of weakly ion- 
ized species to explain the FIP effect, could produce sig- 
nificantly stronger isotope effects. This is concluded from 
analogy with the severe isotope fractionation observed in 
energetic particles related with impulsive flare events. 

Inefficient Coulomb drag is most probably the most effi- 
cient isotopic fractionation process. Under the assumption 
of an isothermal solar atmosphere and neglecting wave ac- 
celeration, one may derive the following expression for the 
maximum fractionation between the two isotopes k and I of 
an clement due to inefficient Coulomb drag [Bodmer, 1996; 
Bodmer and Bochsler, 1998; Bodmer and Bochsler, 2000] 



fki oc 



2Afc-g-i Ak + i 



2Ai-Q-iy Ai + 1 



This evaluates to more than 5% per mass unit for the ex- 
ample of silicon. Thus, for isotopic fractionation, inefficient 
Coulomb drag can be two orders of magnitude larger than 
the FIP effect. The fractionation of isotopes has been tested 
for in the solar wind. Combining measurements of Ne, Mg, 
and Si, Kallenbach et al. [1998b] found a depletion of the 



heavier isotopes in the slow solar wind. However, the iso- 
topic composition of the solar wind is consistent with me- 
teoritic values in high-speed streams [Bochsler et al., 1997, 
Kallenbach et al., 1997c]. This confirms theoretical expec- 
tations that isotopic fractionation should be very weak in 
high-speed streams [Bodmer and Bochsler, 1998; Bodmer 
and Bochsler, 2000]. Since the proton flux is high in the 
high-speed solar wind, Coulomb friction is probably unim- 
portant in these streams, consistent with UVCS measure- 
ments of high outflow velocities of heavy ions in coronal 
holes [e. g. Kohl et al, 1997] which are the source regions of 
the high-speed streams [Krieger et al., 1973]. Since these 
observations show oxygen out-flowing at higher speeds than 
protons [Cranmer et al., 1999], it appears that Coulomb 
drag plays no role in carrying heavy species in the high-speed 
solar wind. If indeed wave-particle interaction succeeds to 
impart the momentum needed to move these species into 
the coronal hole associated solar wind, as is generally be- 
lieved, isotope effects are expected to be weak. This is also 
borne out by the observed constancy of the "'He/'^He iso- 
topic abundance ratio in high-speed streams [Bodmer and 
Bochsler, 1998; Gloeckler and Geiss, 1998]. It is thus to 
be expected that the isotopic composition of the fast solar 
wind most accurately resembles that of the reservoir that it 
is being fed from. Hence it is measurements in high-speed 
streams which should be used to limit the amount of gravi- 
tational settling of the heavy isotopes toward the solar core 
during solar evolution. 

4.5. Selection 

In view of the variations in meteoritic elemental abun- 
dance ratios and in the expected magnitude of solar wind 
fractionation processes, it appears unrealistic to use elemen- 
tal abundance ratios as a test for solar evolution models. 
Nevertheless, for completeness sake, we give a list in Ta- 
ble 2 of those elemental abundance ratios for low-FIP (< 10 
V) elements for which the effect of solar evolution exceeds 
2.5 permill. Of these Ca/Mg is the top candidate from the 
evolution point of view, but, as already mentioned, is not de- 
termined with sufficient accuracy in CI chondrites. If some 
clear-cut criterion for sample selection could be found that 
would narrow down the sample variability, this could be the 
best compositional test of solar evolution models. 

On the other hand, the small variations in the isotopic 
composition of refractory elements (c. g. Mg, Si,. . . ) obvi- 
ous from the already discussed Figmo 6 should allow such a 
test. Current in-situ mass spectrometers are not capable of 
measuring the isotopic composition of the solar wind to the 
precision necessary for this study. Linear time-of-fiight mass 
spectrometers do not resolve isotopes of refractory elements 
and isochronous mass spectrometers have systematic uncer- 
tainties which are about two orders of magnitude too large 
to detect the effect of solar evolution. However, the recently 
launched Genesis mission will collect and return to Earth 
samples of the solar wind in various solar wind regimes in- 
cluding the fast solar wind. These samples will be available 
to laboratory analysis. The Genesis mission has the aim of 
measuring the isotopic composition of oxygen to a relative 
precision of lO^**. The abundance of the rare isotopes of 
oxygen are low (^^O/^'^O ~ 1/2625, ^'^O/^^O ~ 1/499) in 
solar system samples [Anders and Grevesse, 1989; Clayton, 
1993] and in the solar wind [Collier et al., 1998; Wimmer- 
Schweingruber et al., 2001]. Because of the large abundance 
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Table 2. Some of the highest photosphciic abundance ratios as 
predicted by solar models for elements of FIP < lOeV. The ab- 
solute values of the difTcrencc in mass and FIP as well as of the 
relative change in abundance ratios at the solar age. The numeral 
on the leftmost column indicates the position in the list of all 
relative variations of the ratios in inverse numerical order, 
order A/B A(mass) A(fip) A(Xa/Xb) (%) 

(amu) (cV) No mixing Mixing 
~ Ca/Mg 15l8 OS TTg 08^ 

4 Ca/Al 13.3 0.13 1.11 0.78 

5 Ti/Mg 23.6 0.83 1.11 0.78 
7 Fe/Ca 15.8 1.76 1.05 0.74 
9 Ti/Al 21.1 0.83 1.03 0.72 

11 Ca/Si 12.0 2.04 0.98 0.69 

12 Fe/Ti 8.0 1.05 0.97 0.68 
14 Ca/Na 17.1 0.97 0.92 0.65 
16 Ti/Si 19.8 1.33 0.90 0.63 
18 K/Mg 14.8 3.31 0.89 0.62 

20 Cr/Mg 27.8 0.88 0.88 0.61 

21 Mn/Ca 14.9 1.32 0.85 0.59 
25 K/Al 12.3 1.65 0.81 0.57 
29 Mn/Ti 7.1 0.62 0.77 0.54 
31 Fe/K 16.7 3.53 0.76 0.53 
36 Fe/Cr 3.9 1.10 0.74 0.52 
51 Ni/Fe 2.8 0.24 0.58 0.41 



of the rare isotopes of Mg and Si isotopes, this measurement 
should be possible with these samples. This measurement 
precision for O isotopes translates to a similarly high pre- 
cision for Mg and Si isotopes. Because the heavy and less 
abundant isotopes of Mg and Si are not nearly as rare as 
those of O, Genesis should return high-precision and statis- 
tically significant samples of the isotopic composition of Mg 
and Si in the solar wind. 

Mg and Si have other positive aspects to them that make 
them attractive for this study. The are the lightest refractive 
elements with three stable isotopes, moreover, their rare iso- 
topes are still fairly abundant. In addition, as can be seen in 
Figure 3, thoy lie in a local minimum in the predicted pho- 
tospheric abundance AXgurface, i- e. a local maximum in the 
change. Because the isotopic effect is well approximated by 
AXsurface times the relative mass difference of the isotopes, 
the isotopes of these two elements offer the best chance of 
measuring subtle changes in their composition due to solar 
evolution. 

5. Discussion 

Today standard solar models include effects of elemental 
segregation. In this picture, all elements heavier than hy- 
drogen are gradually depleted in the solar photosphere over 
the course of solar evolution. Their current photospheric 
abundances are about 8% lower than in the bulk Sun or in 
the presolar nebula. This depletion is primarily due to dif- 
ferences in the radiative and gravitational forces acting on 
the various elements. Such models compare favorably with 
sound speed profiles of the Sun inferred from helioseismol- 
ogy. Nevertheless, the agreement is far from satisfactory, 
differences in inferred and model sound speeds arc still very 
large compared to the modeled uncertainties on their mag- 
nitudes. 

In this paper we have investigated a more direct ap- 
proach to obtain limits on the amount of element segregation 
throughout solar history. The smallness of the effect makes 
it very difficult to detect. There are two possibilities to check 
such a secular change in photospheric abundances. One can 
measure implanted solar wind in lunar soils of various an- 
tiquity, or one can compare contemporary photospheric or 



solar wind abundances with bulk solar abundances inferred 
from (predominantly) meteoritic data. The first test has 
proven to remain ambiguous [Heber, 2002] . Uncertainties in 
photospheric abundance determinations are on the order of 
or exceed 25% - 30%. These uncertainties arc mainly due to 
badly known atomic properties [Del Zanna et al, 2001], but 
also due to the badly understood physics of the photosphere 
and lower chromosphere [Holweger, 2001; Del Zanna et al., 
2001]. Hence it is not to be expected that photospheric 
abundances will be measured in the near future with suffi- 
cient accuracy to measure the effect of element migration in 
the near future. Therefore, we have investigated the possi- 
bility of detecting or at least putting limits on this effect by 
using in-situ solar wind data. 

The use of in-situ solar wind data immediately restricts 
the set of elements than may be considered. Because of the 
magnitude of the insufficiently understood FIP effect, only 
elements with low FIP (FIP < 10 V) should be considered 
for comparison. However, in this case, no comparison can be 
made with hydrogen, implying that the magnitude of the ef- 
fect of element settling is decreased considerably. It amounts 
to only typically 0.5% for most element combinations. We 
give a list of the best candidate pairs in Table 2. The small- 
ness now results in an additional difficulty, namely that the 
normalizing meteoritic abundances (or element ratios) are 
not known with sufficient accuracy. 

We have shown that only isotopes of refractory elements 
offer any chance to measure the magnitude of the effect of 
elemental segregation in the solar convection zone. With 
Genesis it should be possible to determine the Mg and Si 
isotopes with an accuracy of about 1 permill, good enough 
to detect this effect because in this case, meteoritic data are 
of comparable or better quality. 

Alternatively, high-precision measurements of solar wind 
noble gas isotopic abundances in lunar soils still has the po- 
tential to unravel solar chemical evolution if isotopic and 
elemental fractionation in the implantation process can be 
reliably modeled and if it can be understood in sufficient de- 
tail throughout the complex history of lunar soils and other 
archives for the solar wind. 

The alteration of the photospheric (and hence solar wind) 
isotope composition is especially important for oxygen which 
will be measured at a very high precision (~ 0.01%) with 
Genesis. Here, accurate measurements of isotope ratios of 
refractory elements such as Mg and Si will be important to 
"calibrate" the influence of solar evolution on more volatile 
and chemically reactive elements such as oxygen. Inclusion 
of element migration throughout solar history will be cru- 
cial to the interpretation of the measured values of ^^^0/^*0 
and ^^O/^^O. Because solar isotopic evolution will tend to 
follow the traditional chemical fractionation line slope of 
~ 1/2 while alternative models for O isotopes in the solar 
system involve non-mass-fractionating processes, a precise 
knowledge of the solar effects is of utmost importance for 
the understanding of the history of the solar system. 

While it is doubtful that the required understanding of all 
fractionation processes in the solar wind exists at this time 
for an accurate analysis of isotopic ratios in relation to phys- 
ical processes in the Sun but we believe that the necessary 
knowledge is attainable. 
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